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Recently we proposed an extension to the traffic model of Aw, Rascle and Greenberg. The 
extended traffic model can be written as a hyperbolic system of balance laws and numerically 
reproduces the reverse A shape of the fundamental diagram of traffic flow. In the current work we 
analyze the steady state solutions of the new model and their stability properties. In addition to 
the equilibrium flow curve the trivial steady state solutions form two additional branches in the 
flow-density diagram. We show that the characteristic structure excludes parts of these branches 
resulting in the reverse A shape of the flow-density relation. The upper branch is metastable against 
the formation of synchronized flow for intermediate densities and unstable for high densities, whereas 
the lower branch is unstable for intermediate densities and metastable for high densities. Moreover, 
the model can reproduce the typical speed of the downstream front of wide moving jams. It further 
reproduces a constant outflow from wide moving jams, which is far below the maximum free flow. 
Applying the model to simulate traffic flow at a bottleneck we observe a general pattern with wide 
moving jams traveling through the bottleneck. 

PACS numbers: 89.40.Bb, 05.10.-a, 47.20.Cq 



I. INTRODUCTION 

Modeling vehicular traffic flow using methods from 
continuum fluid dynamics has a long history QII2>OJ0>O1> 
with many contributions mainly from traffic engineers, 
physicists and mathematicians. The existing models can 
be subdivided into first and higher order models ac- 
cording to the highest derivative appearing in the par- 
tial differential equations describing traffic flows 0. |2^| 
First order models, as the classical model of Lighthill, 
Whitham and Richards 7, 8] , approximate the higher or- 
der models by neglecting diffusion terms in the same way 
the Euler equation approximates the Navier-Stokes equa- 
tion. As a consequence, discontinuous data can develop, 
which have to be dealt with by a numerical algorithm. 

In the existing literature on first order systems the role 
of a source term in the partial differential equations de- 
scribing vehicular traffic flow on a road section without 
entries and exits has not been studied in depth. Instead, 
the main focus has been laid on the principal part of the 
equations |9(, be. the collection of terms in the partial 
differential equation containing derivatives of order equal 
to the order of the partial differential equati on JlCl. and 
systems with constant relaxation time |TlL TT3 . fl3t.ll4l. fl5| . 
In 0| , we presented the balanced vehicular traffic model 
(BVT model) , which gen eralizes the model of Aw, Ras- 
cle and Greenberg [{J, lllj by prescribing a more general 
source term subsumed under an effective relaxation coef- 
ficient. Unlike in earlier studies, this effective relaxation 
coefficient depends on both traffic densit y a nd velocity. 
As we showed in numerical simulations in [lq , the model 
can reproduce the observed reverse A shape of the fun- 
damental diagram of traffic flow |l7j . 

In the current work we aim at obtaining a better under- 
standing of the numerical results of the BVT model. In 
particular we study the appearance of the new branches 



in the fundamental diagram in the congested regime, 
which finally form the reverse A shape. Our explana- 
tion of the reverse A shape differs from earlier explana- 
tions, which explained the reverse A as a part of a single 
equilibrium flow branch. In contrast, as we will show, 
the (meta-) stable curve sections of three steady state 
branches form the reverse A in the BVT model. Hence, 
the criticism of the one-dimensionality of steady states 
of fluid-dynamical models [j| does not apply to the BVT 
model. 

In the BVT model traffic flow is described by the fol- 
lowing system of balance laws determining the density 
p = p(t, x) and velocity v = v(t, x) of vehicles 



dp d(pv) 
dt dx 
d(p{v - u(p)) d(pv(v - u(p))) 
dt dx 

(3{p,v)p{u{p) - v). 
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As usual, (t, x) denote the time and space variable. u(p) 
denotes the equilibrium velocity, which fulfills 



u'(p) < for < p < p 7l 



d 2 (pu(p)) 
dp 2 



< for < p < p ri 



(3) 
(4) 



The effective relaxation coefficient (3(p, v) fulfills 

f3(p, v) < for < pi < p < p 2 < p m ,v = u(p), (5) 
/3(p, v) > for < p < pi or p 2 < p < p m ,v = u(p), (6) 

lim p(p,v) >0. (7) 

v^0,u m =u(0) 

Note that in moving observer coordinates, Eqn. J5J re- 
duces to 



dt 



- u(p)) = -P(p,v)(v - u(p)), 



(8) 
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i.e. j3 can be interpreted as decay parameter. As an effec- 
tive parameter, which takes into account the actual relax- 
ation time but also the reaction time, the parameter can 
become negative for intermediate to high densities © 
(see |l6|]). Note that for f3 < 0, the average velocity v 
departs from the equilibrium velocity u{p), i.e. \v — u(p)\ 
increases with time. For (3 > 0, drivers approach the 
equilibrium velocity u(p) with a rate determined by the 
value of (3. 

In |l6j ] we used parameter functions and values which 
describe traffic flow only qualitatively. In the current 
work we use the equilibrium velocity of Newell |l8j ] 



u(p) = u m ( 1 - exp f f 

- V U m \p Pro 



(9) 



with parameter values u m = 160 km/h, A = 
3600 [1/h/lane], p, n — 160 [1/km/lane] and an effective 
relaxation coefficient 



0(M = < 



u—v ' 

da 



(3(p,v) 



if P(p, v)(u(p) - v) - a c > 0, 
ifp(p,v)(u(p)-v)-d c <0, 
else, 



(10) 



0(p,v) = J—(\u{p)-v + ai f±v{p)\+a 2 f\v{p)) (If i 

V™ 



and 



Av(p) = tanh (a 3 — ) (u(p) + cp m (^ 
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with parameters a c — 2 m/s 2 , d c — —5 m/s 2 , T = 0.1 s, 
«! = -0.2, a 2 = -0.8, a 3 = 7 and c = -14 km/h. The 
density values, which determine the sign of f3 according 
to Eqs. JSJ-JfjJ), are p\ — 19.09 [1/km/lane] and p 2 = Pm- 
Throughout this work, we model two-lane sections of a 
highway without entries and exits with the above param- 
eter set. In comparison to 0] theses parameters describe 
traffic flow more realistically, although we have not used 
experimental traffic data to determine them for a specific 
highway section. Note that the general conditions JSJ)- 
(|7|) are sufficient to obtain multivalued fundamental dia- 
grams. However, the quantitative details depend on the 
precise choice of (3. As our simulations show, the analyt- 
ically derived properties of the model of Aw, Rascle and 
Greenberg, that the velocity does not become negative 
and collisions do not occur, carry over to our system. For 
the numerical simulations of the model equations Q-© 
we used a high-resolution shock-capturing scheme with 
an approximate Riemann solver. We chose a spatial res- 
olution of 20 m and dynamically adapted the temporal 
resolution to half the value obtained from the Courant 
condition. The numerical method is described in detail 
in [l|. 

In order to obtain a deeper insight into the structure 
of the BVT model we study the smooth steady state so- 
lutions in Sec. ^ Our numerical simulations produce 



more general solutions approximating steady states so- 
lutions, which will be discussed in Sec. 11111 In Sec. II VI 
we assess the stability properties of the steady state so- 
lutions. With theses results we classify the traffic states 
of the BVT model according to the three traffic phases 
of Kerner [j| in Sec. [V] and apply the BVT model to sim- 
ulate traffic flow at a bottleneck in Sec. lVIl We conclude 
the paper in Sec. IVIII 



II. SMOOTH STEADY STATE SOLUTIONS OF 
THE BVT MODEL 

For smooth solutions the balance equations describing 
traffic flow 0-10) can be rewritten as 



dp dp dv 
dt dx ^ dx 
dv , dv 

m + {v + pu {p)) d-x 



0, (13) 
(3(p,v)(u(p)-v). (14) 



In the following, we study the smooth steady state solu- 
tions of the BVT model. In comparison to the study by 
Lee, Lee and Kim [l9l | , our analysis - being performed on 
a first order system - is considerably simpler. Note, how- 
ever, that due to the possibility of dealing with discon- 
tinuous solutions, we can in principle have more general 
solutions in the balanced system CQ)-I0 (see also [l( and 
references therein). 

Let us repeat, that for a steady state solution, there is 
a coordinate system (i, z) and a constant velocity w 



x 
t 



wt, 



such that 



dt 

d -l = 0. 

dt 



(15) 
(16) 



(17) 
(18) 



It follows from the continuity equation 113|) that for all 
steady state solutions there is a constant q such that 



pv = q + pw. 



(19) 



Hence steady state solutions are restricted to straight 
lines in the fundamental diagram of traffic flow. More- 
over, the minimum and maximum speed of information 
propagation in system (Q-© limit the physically admis- 
sible steady state solutions, i.e. the velocity w, as 



Ai = v + pv! (p) < w < A2 = v. 



(20) 



Let us assume that v — w > (else q = 0) . Then we can 
solve Eqn. (|19|) for the density 



(21) 



and substitute it into Eqn. I|14|l . yielding the ordinary 
differential equation for steady state solutions in the BVT 
model 

(Ax - «) ^ = - «). (22) 

\ J dz v — w \ v — w ) 



A. Trivial steady state solutions 

We first study the trivial (constant) solutions of this 
ODE, i.e. solutions fulfilling 4| = 0. The solutions are: 

• The equilibrium velocity curve: 

v = v e = u{p). (23) 

• The jam line (compare to the line J of Kerner pj): 

v = v ] = u{p) + (ai + a 2 )Av(p) (24) 
for pi < p < p 2 . 



Note that pv < pu(p). 
• The high-flow branch: 

v = v h = u(p) + (a 



i - a 2 )Av(p) 
for pi < p < p 2 . 



(25) 
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FIG. 1: (Color online) The trivial steady state solutions (i.e. 
solutions with 4^ = 0) of the BVT model. These solutions are 
the equilibrium solution v = u{p) (solid black curve) and the 
two branches fulfilling fl(p,v) — (dashed curves). Moreover, 
we present the results of simulation runs of perturbed equi- 
librium data. The simulated data points in the fundamental 
diagram are closely related to sections of trivial steady state 
solutions fulfilling (3{p, v) — 0. 



Note that pv > pu(p) in this case. 

We summarize the trivial steady state solutions in Fig.^ 
We further show in this figure the results of simulation 
runs for constant initial data in equilibrium p = po, 
po = 1,2,. ..,159 [1/km/lane], v = u(po) on a 7 km long 
stretch of a highway with periodic boundary conditions, 
prescribing a small amplitude perturbation of the den- 
sity bp = sin(7rx) initially located between 2 and 3 km 
on top. The data points were extracted from the simula- 
tions at 5 equidistantly distributed virtual detectors after 
an evolution time of 10 h, without applying a temporal 
aggregation. As one can see from the plot many data 
points of the numerical solutions are closely related to 
branch sections of trivial steady state solutions, in par- 
ticular for the jam line. 

B. Non-trivial steady state solutions 

Let us compare the equation of steady states (1221) in 
an alog y to the dynamics in classical mechanics presented 
in [jjj to the following equation 

(Xi- W ) d £ = -jLu( V , W , q) , (26) 

with a potential energy U. We find that the potential 
U(v, ...) has a functional form which is camelback- shaped 



for a wide range of constant values w and q as in |19j . 
However, unlike in |19|. we do not obtain an accelera- 
tion term in our first order system. Note that according 
to fflty Ai - w < 0. 

In the following we restrict the discussion to the phys- 
ically admissible smooth solutions (in particular we do 
not consider solutions with infinite gradient). In the limit 
z — * ±oo all maximally extended steady state solutions 
approach one of the following curve sections: 

A: free equilibrium flow: v = u(p) and p < pi, 

B: unstable equilibrium flow: v = u(p) and p\ < p < p 2 , 

C: jam line: see Eqn. I|24l) . 

D: high- flow branch: see Eqn. I|25(l . 

We classify the non-trivial maximally extended steady so- 
lutions according to the behavior in the limit z —* ±oo, 
using the letters of the corresponding branches. Due 
to the characteristic structure of the BVT model (see 
Fig- El 7 on ly five solution classes can appear. These are 
the classes AD, BC, BD, CC and DD j2g. Schematically, 
all these steady state solutions have the form indicated 
in Fig. i.e. solutions lying in the regions II and IV of 
Fig. fulfill the condition > 0, solutions in regions I 
and III obey 4^ < respectively and lim z ^± oc ^ = 0. 
Comparin g ou r classification to the classification of Lee, 
Lee, Kim |l9j for cases, where the parameters q and w 
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FIG. 2: (Color online) Characteristic structure and steady 
state solutions of the BVT model. As solid lines we plot 
the trivial steady state solutions bordering the regions I - 
IV. On top (dashed curves) we plot the characteristic curves 
with slope Ai and A2 respectively. As the speed of steady 
state solutions is limited by the characteristic speeds Ai and 
A2, physically admissible steady state solutions lie inside the 
characteristic cones spanned by these two speeds at every 
point (p,pv) in the flow-density diagram. 



FIG. 3: Sketch of the non-trivial steady state solutions. The 
steady state solutions fulfill > (or 4^ < respectively) 

dv 
dz 



lead to a camelback-shaped profile of the potential U, the 
solutions correspond to minimum-saddle solutions (with- 
out oscillations). The possible appearance of the different 
solution classes in the flow-density diagram is shown in 

Fig. m 



III. QUASI STEADY STATE SOLUTIONS 

In principal, some of the non-trivial steady state so- 
lutions described before can be glued together to form 
discontinuous, periodic steady state solutions [Til fl5| . 
Steady state solutions can be linked by a shock wave, if 
the quantities p_, u_ left to the interface and the corre- 
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FIG. 4: (Color online) Left Panel: Regions covered by the 
smooth maximally extended non-trivial steady state solutions 
linking different branches, which are consistent with the char- 
acteristic structure of the BVT model. The solutions link 
the free equilibrium flow solution with the high-flow branch 
(class AD), the unstable equilibrium solution with the jam 
line (class BC) and the unstable equilibrium solution with 
the high-flow branch (class BD). 

Right Panel: Regions covered by the smooth maximally 
extended non-trivial steady state solutions linking identical 
branches, which are consistent with the characteristic struc- 
ture of the BVT model. For class CC the non-trivial steady 
states link data points lying on the jam line, for the class DD 
they link data points of the high-flow branch. 



sponding quantities p+, v+ rig ht of the interface satisfy 



the following conditions 



Sng 



P+ > P-, 

p+(v + - u(p+)) = p-(v- - u(p-)). 
With the equation of steady state solutions 

P±v± =q + p±w, 
we obtain from Eqn. (|28|l 

w{p+ - P-) = P+u{p+) - p-u(p-). 



(27) 



(28) 



(29) 



(30) 



The monotonicity of the velocity of steady state solu- 
tions, which follows from Eqn. I|22(l (see also Fig. 03) , re- 
lates to the monotonicity of the density according to 



dp 

dz 



dv 



w dz 



(31) 



Therefore, condition (|2Tjl restricts periodic steady state 
solutions linked by shock waves to regions II and IV of 
Fig. El 

In the following we focus on steady state solutions of 
type CC in region II. In particular, steady state solutions 
lying inside region II have to approach the jam line at 
velocity values for which the condition 



(pvi(p))">0 



(32) 
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is fulfilled. For the chosen parameter values, this implies 
p_ > 36.51 [1 /km/lane]. A necessary requirement to 
fulfill Eqn. (j20J is the condition 

{pu{p))'\ p=p _ > (Mp)W- (33) 

which results in our particular case to p- < 61.57 
[1/km/lane]. Therefore, for the chosen parameter values 
there is only a very small parameter range where periodic 
steady state solutions of type CC linked by shock waves 
are possible. However, in our numerical simulation, we 
find quasi steady state solutions in a much larger region of 
parameter space (see Fig.^J. To obtain a more detailed 




FIG. 5: Quasi steady state solution at the jam line for a nu- 
merical simulation of perturbed equilibrium data with initial 
density p — 100 [1/km/lane]. Between a simulation time of 8 
and 10 h the evolution is quasi-stationary with a propagation 
speed of about w ~ — 14 km/h. 

picture of these quasi stationary solutions, we plot the 
simulation results of unstable equilibrium data in Fig. 
focusing on the evolution between 8 and 10 hours simula- 
tion time. We started the simulation of Fig. El with con- 
stant initial data p = po = 100 [1/km/lanc], v = u(po) 
and a sinusoidal density perturbation on top. After a 
long evolution time we obtain a quasi steady state close 
to the jam line, propagating with a velocity w ~ — 14 
km/h upstream. This quasi steady state solution is not a 
true steady state solution, the amplitude slowly decreases 
with time, as it can be scarcely noted in the figure. The 
solution consists of different branches close to the triv- 
ial steady state solutions lying in regions I /II of Fig. 
Thus, for quasi steady state solutions the quantities q 
and w of steady state solutions (see Eqn. tylty ) are only 
approximately constant. 

Similar quasi-steady state solutions also exist for a sec- 
tion of the high- flow branch (see Sec. IIVI below). They 
consist of approximate steady state solutions lying in re- 
gions III/IV of Fig. |21 These quasi steady state solutions 
travel downstream with a velocity approximately given 
by the tangent of the high-flow branch. 



IV. STABILITY ANALYSIS OF THE STEADY 
STATE SOLUTIONS 

Due to the importance of steady state solutions at the 
jam line and the high-flow branch (see Fig. ^) we focus 
the stability analysis on the trivial steady state solutions 
of the BVT model. In principal the code can be applied 
to study the stability of the non-trivial solutions as well, 
although more accurate results may be obtained using 
specialized methods to this aim (see e.g. [2(J ). 

A. Linear stability analysis 

As presented in 0] the equilibrium flow curve pv = 
pu{p) is linearly stable for p < p\ and linearly unstable 
in the intermediate to high density regime p± < p < 
P2- Here we extend the linear stability analysis to all 
steady state solutions obtained from setting (3(p, v) = 
0, i.e. steady state solutions Eqs. $M and We 
denote the corresponding constant states (po,v^ h (po))- 
Plugging the ansatz 

p = po + pexp(ilx + w(l)t), (34) 
v = v j/h (p )+vexp(ilx + ui(l)t) (35) 

into the evolution equations (|T3)) and (fHI) we obtain 
as characteristic equations for the existence of solutions 
(p,v) ^ (0,0): 

(u + ilv j/h ) 2 + (w + ilv j,h )(ilu' Po + (v i/h -u)^-) 

ov 

-p i&v^ h -u)=0. (36) 
dp 

Solving the last equation for u> we can distinguish be- 
tween the linearly stable and unstable regime of the 
trivial steady-state solutions, i.e. curve sections with 
Re(w) < for arbitrary / and curve sections for which 
Re(w) > for some I respectively. We find that the 
jam line is linearly unstable for pi < p < pj = 39.73 
[1/km/lane] and linearly stable for pj < p < p2, whereas 
the high-flow branch is linearly stable for densities p\ < 
P < Ph — 39.73 [1/km/lane] and linearly unstable for 
p h < p < p 2 . 

In the following we will give a more intuitive explana- 
tion for the above results. For trivial steady state solu- 
tions with exactly constant density and velocity profile, 
the characteristic structure H'20|) does not give any restric- 
tions. However, quasi steady state solutions with non- 
constant density (velocity) on these trivial steady state 
branches travel with a velocity w which corresponds to 
the derivative of the flow density curve 

*n«Bp.. (37) 
dp 

Hence, according to condition l|20|) the local character- 
istic cone has to enclose the corresponding steady state 
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branch spanned by the characteristic speeds Ai and A2, 
i.e. Ai < w < A2. For our parameter values the 
characteristic condition restricts possible stable (quasi) 
steady state solutions at the jam line to solutions fulfilling 
p > jjj = 39.73 [1 /km/lane] and those at the high-flow 
branch to solutions fulfilling p < ph = 39.73 [1/km/lanc] 
(see Fig. |2J). Note, that the two densities need not agree 
for a general parameterization. For our parameteriza- 
tion, it follows from the appearance of the term Av(p) in 
both, Eqs. and 

Next we study the nonlinear stability properties in sim- 
ulations of the full system. 
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FIG. 6: (Color online) Left panel: Stability properties of 
the high-flow branch. We study the stability properties by 
prescribing constant steady state solutions p = po, v = 
v h (po) and on top a sinusoidal velocity perturbation Sv = 
t>ampl sin(vT3;) for 2 < x < 3 km. In the figure, we plot 
the maximum value of v h (po) + Wampi, Wampi < 0, for which 
the initial data is unstable against the formation of synchro- 
nized flow, for different evolution times. For densities p < 40 
[1/km/lane], the high-flow branch is metastable, it becomes 
unstable against the formation of synchronized flow for suffi- 
ciently small values « anlp i. For densities p > 40 [1/km/lane] 
the high-flow branch is unstable against the formation of syn- 
chronized traffic flow. 

Right panel: Analysis of the stability properties of the jam 
line using constant steady state solutions p — po, v — 
v 3 {po) and on top a sinusoidal velocity perturbation Sv — 
Wampi sin(vTx) for 2 < x < 3 km. In the figure, we plot the min- 
imum value of v 3 (po) + «am P i, Wampi > 0, for which the initial 
data is unstable against the formation of synchronized flow, 
for different evolution times. For densities p > 40 [1/km/lane] 
the jam line is metastable, it becomes unstable only for suf- 
ficiently large values « amp i. For densities p < 40 [1/km/lane] 
the jam line is unstable against the formation of synchronized 
traffic flow. 



B. Numerical results 

We first analyze the stability properties of the high- 
flow branch of steady state solutions, Eqn. To this 
aim, we use constant steady state initial data p — po , v — 
v (po) with a sinusoidal perturbation Sv — v amp \s'm(irx) 
for 2 < x < 3 km, prescribing periodic boundary con- 
ditions on a 7 km long highway. In order to decide 



whether synchronized flow appears during the numeri- 
cal evolution, we use the criteria v(p) < u(p) or v(p) < 

V h (p) - |Uampl|< 

As our analysis shows, the stability properties of the 
high-flow branch depend on the particular perturbation. 
For the density regime pi < p < Ph, the high- flow branch 
is metastable against the formation of synchronized flow, 
i.e. for small amplitude perturbations, no synchronized 
flow appears, whereas for larger velocity perturbations 
with negative amplitude w am pi, synchronized flow ap- 
pears. For the density regime ph < p < pi-, the high- flow 
branch is unstable against the formation of synchronized 
flow. We summarize the corresponding results in the left 
panel of Fig. HJ1 

Second, we study the stability properties of the jam 
line, using constant steady state initial data p = po, 
v = v J (po) again with a sinusoidal perturbation Sv — 
Uampl sin(yTiz;) for 2 < x < 3 km. In this case, we use the 
criteria v(p) > u(p) or v(p) > v J (p) + \v amp \\ to identify 
synchronized flow. 

For densities pj < p < p2 the jam line is metastable 
against the formation of synchronized flow, i.e. for small 
amplitude perturbations, the jam line is stable, for larger 
amplitude perturbations (with positive amplitude v amp \), 
synchronized flow appears. In contrast, for densities 
pi < p < pj [1/km] the jam line is unstable against the 
formation of synchronized flow. We show the results in 
the right panel of Fig. 

We summarize the results of the stability analysis in 
Fig. The figure shows the observed gap between free 
and congested flow in the fundamental diagram of traffic 
flow separating the (meta-) stable branch sections. 



V. IDENTIFICATION OF KERNER'S THREE 
PHASES OF TRAFFIC FLOW 

Kerner 0, 121L l22l |23| classifies traffic flow into three 
phases: free flow, synchronized flow and wide moving 
jams. In this section we try to relate the traffic states of 
the BVT model to Kerner's three phases. We will first 
summarize our classification, before we will discuss the 
motivation. 

• free flow: Steady state solutions at the equilibrium 
velocity curve u = u(p) for the density regime < 
P < Pi (free equilibrium flow) and (quasi) steady 
state solutions close to the high-flow branch in the 
metastable regime pi < p < ph make up the free 
flow state. 

• wide moving jams: Spatially extended (quasi) 
steady state solutions at the jam line in the 
metastable regime pj < p < p2 make up wide mov- 
ing jams. 

• synchronized flow: All other congested traffic states 
including the non-trivial steady state solutions of 
type BC and BD form synchronized traffic flow. 
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FIG. 7: (Color online) Results of the stability analysis of the 
trivial steady state solutions 4^=0. Curve sections repre- 
sented as solid black line correspond to the linearly stable 
steady state solutions. Curve sections represented as dotted 
line and dashed line correspond to unstable and metastable 
steady state solutions respectively. 



A. Free flow 

For small densities (0 < p < p\) free flow is stable in 
the BVT model. Moreover, the model can reproduce the 
metastability of free flow against the formation of syn- 
chronized flow (see Fig. , which is observed for traffic 
states at the high- flow branch p\ < p < ph- In the model 
instabilities only appear for velocity perturbations with 
negative amplitude. 



B. Wide moving jams 

As our results show wide moving jams as defined above 
are stable against small amplitude perturbations. For the 
chosen parameter values the propagation speed of wide 
moving jams lies between —16 < w < —14 [km/h] for the 
density region p t < p < p2, i.e. it is nearly constant and 
reproduces the observed value. 

We further analyzed the outflow from wide moving 
jams (see Fig. |SJ. We find for wide moving jams in the 
density region p t < p < 125 [1 /km/lane] a constant out- 
flow of about / ou t = 1914 [1/h/lane]. This can be seen 
in the lower right panel of Fig. [HI where we plot the out- 
flow from wide moving jams p — Pj,v = v 3 (pj) initially 
located between 2 and 3 km and surrounded by a region 
of free flow p = pt,v — u(pf). 
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FIG. 8: (Color online) Upper two panels: Simulation of the 
outflow from a wide moving jam. Plotted are the evolution of 
the density and the velocity as a function of space and time. 
In the initial data we prescribe a wide moving jam with den- 
sity pj — 100 [1/km/lane] between 2 and 3 km, and free flow 
with density pf — 5 [1/km/lane] elsewhere. During the evolu- 
tion, the wide moving jam narrows down and finally dissolves. 
Lower left panel: Flow-density diagram for the above simu- 
lation results after t = 0.02 h. In addition we plot the curve 
representing the equilibrium flow as well as the jam line and 
the high-flow branch. We determine the outflow of the wide 
moving jam at that point for which the velocity of the outflow 
differs from the equilibrium velocity by less than 1 %. 
Lower right panel: Outflow from wide moving jams, varying 
the jam density pj of the wide moving jams between 2 and 3 
km. The three different curves correspond to different values 
of free flow density p/ in the region between and 2 km and 
3 and 7 km. As one can see from the plot, the outflow from 
the wide moving jams only varies within a very small range 
of flow values and it is largely independent of the density of 
free flow. The typical outflow from wide moving jams for 
the chosen parameter values is / out ~ 1914 [1/h/lane], which 
is far below the maximum of metastable free flow / = 2487 
[1/h/lane]. 



C. Synchronized flow 

Synchronized flow as defined above covers a wide re- 
gion in the fundamental diagram. This can be already 
seen in the left panel of Fig. 0] for steady state solu- 
tions of type BC and BD. We exemplarily show the 
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formation of synchronized flow from free flow of den- 
sity p = 30 [1 /km/lane]. As a nucleus for the emer- 
gence of synchronized flow, we use a velocity perturba- 
tion Sv = — 7sin(7ra) located between 2 an 3 km on a 
highway with periodic boundary conditions. The evolu- 
tion of these initial data leads to states which are widely 
scattered in the fundamental diagram, as it can be seen 
in Fig. [5] We remark that in an earlier work [24| the 
scattering was reproduced by modeling different vehicle 
types. Here, in contrast, the scattering already follows 
from the traffic dynamics without distinguishing between 
different vehicle types. Including different vehicle types 
into the BVT model, which is beyond the scope of the 
current work, would be expected to further widen the 
scattering. 
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longitudinal extension of the highway of 7 km with homo- 
geneous initial free flow of density p — 37.5 [1/km/lane]. 
We model the bottleneck simply by a velocity modifica- 
tion (velocity drop for free flow) between 5 and 6 km, 
setting the velocity to a modified value v mod after each 
numerical evolution step, 

v mod = v+ (u(p) - v - 0.1 km/h)| sin(7ra)| (38) 

We show the numerical evolution in Fig. Despite 
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FIG. 10: (Color online) Formation and propagation of wide 
moving jams. The upper panel shows the evolution of the 
density, whereas the lower panel shows the evolution of the 
velocity. At the bottleneck located between 5 and 6 km, syn- 
chronized flow forms, which finally leads to a wide moving 
jam. This wide moving jam moves with a velocity of about 
-15 km/h (i.e. upstream) and swallows moving jams during 
this propagation. It further travels through the bottleneck. 
See the text for a detailed description. 



FIG. 9: (Color online) Formation of synchronized flow from 
metastable free flow of density p — 30 [1/km/lane] with an 
initial velocity perturbation 5v = — 7sin(7ra;) located between 
2 an 3 km. Due to the velocity perturbation, the free flow 
state breaks down, leading to a complicated pattern of syn- 
chronized flow which covers a wide region of states in the 
fundamental diagram. Finally, moving jams form, which can 
in turn lead to free flow of lower density, thus reproducing the 
hysteresis effect observed in traffic dynamics. The plot shows 
all data points corresponding to the constant time slices at 
t = iAt, where At = 0.1 h, i = 0, .., 50. 



VI. TRAFFIC FLOW AT BOTTLENECKS 

In this section we study the behavior of traffic flow in 
the BVT model at a highway bottleneck. We focus the 
discussion on two simulation runs of a two-lane highway 
with periodic boundary conditions |27) |. Again, we use a 



the simplicity of the initial setup the numerical evolution 
shows a very complicated dynamics. As we will discuss 
below, we observe the formation of synchronized flow and 
wide moving jams. 

As a consequence of the bottleneck, the initial veloc- 
ity drops to smaller values in the bottleneck region (dark 
blue regions between 5 and 6 km in the velocity plot), but 
also further upstream (dark blue regions between 4 and 
5 km at about O.lh). Both regions correspond to syn- 
chronized flow. The first synchronized flow region stays 
fixed at the bottleneck, however the upstream front can 
oscillate in time (e.g. between 0.8 h and 0.9 h). The 
second region of synchronized flow travels upstream. It 
takes some time until an accentuated wide moving jam 
with velocities close to zero forms. This wide moving 
jam travels further upstream and reenters the numeri- 
cal domain at 7 km after t ss 0.4 h due to the periodic 
boundary conditions used in the numerical simulation. 
When reaching the bottleneck, it simply travels through 
the first synchronized region, thus becoming a foreign 
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wide moving jam. Note, that the velocity of the down- 
stream front of this wide moving jam is nearly constant 
and has a value of about -15 km/h. 

Between the wide moving jams we observe regions of 
low density and high velocity, which correspond to free 
flow (see e.g. the region at x = 1 km for t = 0.7 h), and 
smaller moving jams. As one can see from the plot there 
are several regions (pinch regions) where these additional 
moving jams form (see e.g. the region between and 3 
km for t R* 0.15 h or at about 1.5 km for t « 0.45 h). For 
these moving jams, the downstream front is in general 
not as robust as for the wide moving jam described above. 
At x = 3.5 km for t = 0.55 h, we observe the merging 
of two moving jams, which are finally swallowed by the 
wide moving jam at about x = 2 km for t = 0.75 h. We 
also observe an example for the catch effect of a narrow 
moving jam, see the region at x = 5.5 km for t w 0.15 h. 
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FIG. 11: (Color online) Formation and propagation of wide 
moving jams. Fhe simulation setup is identical to that of 
Fig. 1101 except that after a time of 0.5 h the bottleneck is 
not effective any longer. As a consequence, the synchronized 
flow region pinned to the bottleneck region disappears, the 
corresponding wide moving jam becomes wider at late times. 



Figure ^2 shows the simulation results for the same 
initial setup, except that the bottleneck Eqn. (|38|l is only 
effective for simulation times t < 0.5 h. As a consequence 
the synchronized flow region pinned to the bottleneck 
disappears at later times. The wide moving jam trav- 
eling through the former bottleneck region expands and 
reaches velocities close to zero inside the jam. This can 
be also seen in Fig. ^] where we plot the time series of 
the flow rate and the velocity for a detector located at 
x = km. 



FIG. 12: Time series of the flow rate (upper panel) and the 
velocity (lower panel) as measured by a detector at location 
x — km for the simulation results of Figs. I10| and llll Inside 
the wide moving jam, small values of the traffic flow and in 
particular the velocity are reached. 



VII. DISCUSSION 

The BVT model is a macroscopic, deterministic model, 
which describes vehicular traffic flow using standard 
methods from continuum fluid dynamics. It uses an equi- 
librium flow-density curve. In contrast to earlier models, 
however, the parameter range of the effective relaxation 
coefficient is extended to negative values. As a conse- 
quence, the equilibrium flow curve does not describe traf- 
fic states in the congested regime directly, but still deter- 
mines the characteristic structure of the model. An ad- 
ditional consequence of the negative effective relaxation 
coefficient is the appearance of (two) additional branches 
of trivial steady states. The characteristic structure, i.e. 
the finiteness of propagation speeds, restricts the stabil- 
ity of these steady state solutions. The high- flow branch 
is metastable against the formation of synchronized flow 
for intermediate densities and unstable for high densi- 
ties. We interpret the metastable section of high-flow 
branch as metastable free flow. Stable free flow, in con- 
trast, corresponds to the stable equilibrium flow. The 
low flow branch in the congested regime (i.e. the jam 
line) is unstable against the formation of synchronized 
flow for intermediate densities and metastable for high 
densities. We interpret spatially extended solutions at 
the metastable branch of the jam line as wide moving 
jams. We further identify the unstable sections of the 
high- flow branch and the jam line, as well as the addi- 
tional (steady state) solutions in the congested regime, 
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which can lead to very complicated oscillatory patterns, 
as synchronized flow. Thus, synchronized flow covers a 
wide region of congested states in the fundamental di- 
agram, without distinguishing between different vehicle 
types and driver characteristics in the model. 

There are some additional results supporting the BVT 
model. In particular, the model ensures that wide moving 
jams do not form spontaneously from free flow. When the 
velocity drops below the critical value in free flow, the 
velocity is driven to even smaller values at that location 
which results in a strong gradient in the velocity. It is 
only after complicated oscillations have occurred, which 
lead to a rearrangement of the density and the velocity, 
that an extended steady state solution close to the jam 
line, i.e. a wide moving jam, appears. 

We further can reproduce the characteristic properties 



of real wide moving jams with our model. For the chosen 
parameter values, the downstream front of wide mov- 
ing jams travels upstream with a nearly constant veloc- 
ity of about 15 km/h. Moreover, the outflow from wide 
moving jams is largely independent of the characteristics 
of the wide moving jam. We obtain a typical outflow 
of 1914 vehicles/h/lane. Furthermore, we showed, that 
wide moving jams travel through bottlenecks, whereas 
smaller moving jams can be caught by a bottleneck. 



Acknowledgments 

We thank James Greenberg for encouraging comments. 



[1] R. Kiihne and P. Michalopoulos, in Transportation Re- 
search Board special report 165, (5)1 (1997, available 
online at http://www.tfhrc.gov/its/tft/tft.html. 

[2] D. Helbing, Verkehrsdynamik (Springer, Berlin, 1997). 

[3] D. Helbing, A. Hennecke, V. Shvetsov, and M. Treiber, 
Transportation Research Part B 35, 183 (2001). 

[4] D. Helbing, Reviews of Modern Physics 73, 1067 (2001). 

[5] B. Kerner, The Physics of Traffic (Springer, Berlin, 
2004). 

[6] F. John, Partial Differential Equations (Springer, New 
York, 1982). 

[7] M. Lighthill and G. Whitham, Proceedings of the Royal 

Society A 229, 317 (1955). 
[8] P. Richards, Operations Research 4, 42 (1956). 
[9] A. Aw and M. Rascle, SIAM Journal on Applied Math- 
ematics 60, 916 (2000). 
[10] P. DuChateau and D. Zachmann, Applied Partial Differ- 
ential Equations (Dover, Mineola, 2002). 
[11] J. Greenberg, SIAM Journal on Applied Mathematics 62, 
729 (2001). 

[12] M. Rascle, Mathematical and Computer Modelling 35, 
581 (2002). 

[13] W. Jin and H. Zhang, Proceedings of the TRB Annual 

Meeting, Washington DC (2003). 
[14] J. Greenberg, A. Klar, and M. Rascle, SIAM Journal on 

Applied Mathematics 63, 818 (2003). 
[15] J. Greenberg, SIAM Journal on Applied Mathematics 64, 

1175 (2004). 



[16] F. Siebel and W. Mauser, SIAM Journal on Applied 

Mathematics 66, 1150 (2006). 
[17] M. Koshi, M. Iwasaki, and I. Ohkura, in Proceedings 

of the 8th International Symposium on Transportation 

and Traffic Theory (V. Hurdle, E. Hauer, G. Stuart ed., 

1983), pp. 403-424. 
[18] G. Newell, Operations Research 9, 209 (1961). 
[19] H. Lee, H.-W. Lee, and D. Kim, Phys. Rev. E 69, 016118 

(2004). 

[20] F. Bouchut, Nonlinear Stability of Finite Volume Meth- 
ods for Hyperbolic Conservation Laws and Weil-Balanced 
Schemes for Sources (Birkhauser, Basel, 2004). 

[21] B. Kerner and H. Rehborn, Phys. Rev. E 53, R4275 
(1996). 

[22] B. Kerner, Physical Review Letters 81, 3797 (1998). 
[23] B. Kerner, Phys. Rev. E 65, 046138 (2002). 
[24] M. Treiber and D. Helbing, J. Phys. A 32, L17 (1999). 
[25] The system studied in this paper is a first order system 

consisting of two equations, see Eqs. Q-Q. 
[26] For the parameters used in [T(| . for which p2 < p m , we 

would have an additional branch section (E) of stopped 

equilibrium flow v = u(p) for pi < p < p m and a sixth 

solution class DE. 
[27] Using periodic boundary conditions enables us to study 

the propagation of moving jams through a bottleneck 

modeling only a single bottleneck. 



